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Abstract: We show how the effects of azimuthal optical aberrations 
on singular light beams can result in an intensity modulation in the beam 
waist or focal point spread function (PSF) that is directly proportional to 
the amplitude of the applied phase aberration. The resulting distortions are 
enough to significantly degrade their utility even in well corrected optical 
systems. However we show that pattern of these intensity modulations is 
related to the azimuthal order of the applied aberration and we suggest how 
this can be used to measure those aberrations. We demonstrate a closed 
loop system using a liquid crystal spatial light modulator as a programable 
diffractive optical element to both generate the beam and correct for the 
sensed aberrations based on feed back from a CCD detected intensity image 
of the focal point spread function. 
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1. Introduction 

Singular light beams contain screw phase dislocations across the wave-front often characterised 
by their topological charge, the multiple of 2/r in phase that the wave traverses along a circular 
path around the dislocation. Such beams, sometimes also referred to as optical vortices or or- 
bital angular momentum (OAM) states dEL have attracted interest in a wide range of fields 
including quantum communication cryptography (H El, optical trapping [6l and far field 
super-resolving optical microscopy [8) . The properties of these beams that make them at- 
tractive for such applications include the sharp zero of intensity that occurs at the dislocation, 
the circularly symmetric donut intensity distribution they result in when focussed and the fact 
that modes with different topological charge can be considered as distinct quantum states of 
light. Unfortunately, when optical aberrations are superimposed on these beams, for instance 
by passing through the atmosphere or simply due to optical misalignments, these properties can 
be severely degraded. For instance, it has been shown that aberrations due to atmospheric 
turbulence results in severe mixing of OAM states, sufficient to destroy the information in a 
quantum encoded data stream. 

In this paper we concentrate on applications where the singular beam is focused with a lens. 
In the ideal case the resulting focal point spread function consists of a bright ring of intensity 
surrounding a sharp zero of intensity. However, we show that this ideal case is quickly degraded 
even in what is normally considered to be well corrected optics. Indeed it is rare to see exper- 
imental results presenting focussed singular beams that are not imperfect or distorted in some 
way. We show analytically that the focal point spread functions of singular beams have off-axis 
linear sensitivity to both odd and even azimuthal aberrations while those of non- singular beams 
have off-axis linear sensitivity only to odd azimuthal aberrations. We utilize this property of the 
singular beam to detect the corresponding aberrations using a segmented detector. We then im- 
plement this technique in a closed-loop system based on feedback from a CCD camera acting 
as the detector, using a liquid crystal spatial light modulator as the beam shaping device. 

2. The sensitivity of singular and normal beams to azimuthal aberrations 

Fig Q] illustrates the effect of aberrations on the point spread function of a simple lens with its 
circular pupil filled with a constant intensity singular beam of varying topological charge, /. 
The intensity point spread function is evaluated in the paraxial scalar limit by calculating the 
modulus squared of the Fourier transform of the pupil function for various applied phase aber- 
rations corresponding to the lowest order Zernike aberration modes. Zernike polynomials Z[ are 
often used to describe optical aberrations over a circular aperture, ifTTIl . and can be expressed 
as a the product of a radially varying component R\(f) and an azimuthally varying component 
Ti(6). In fig Q] an amount of aberration for each Zernike mode has been applied to result in 
a root mean square phase aberration of Qrms = 0.447 radians, sufficient to reduce the Strehl 
ratio of an otherwise perfect lens to 0.8 - the minimum figure considered for a well corrected 
svstem |[T()ll . The resulting PSFs for a planar beam are also shown and clearly exhibit limited 
distortions compared to those visible on the singular beams. For instance astigmatism distorts 
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Fig. 1. Intensity patterns in the Fourier plane for a normal beam (/ = 0) and singular 
beams (I = 1,2,3) aberrated with -0.447Z; (left column) and +0.447Z; (right column) 
of Z i=5 ^ 7 ^ 9 ^ 0M ^ and 15 (from top to bottom). 



the donut intensity pattern into a two lobed pattern. Indeed, astigmatism can be used to con- 
vert between Hermite Gaussian modes and singular Laguerre Gaussian modes using cylindrical 
optics ifHltH. 

The sensitivity of the intensity distribution in a PSF to aberrations can be considered in 
terms of how that intensity varies as the amplitude of the applied aberration is increased. For 
instance, an unaberrated planar beam produces an absolute peak of intensity on-axis that can 
only decrease as any aberration is applied. To maintain an analytic form this means that the 
intensity varies only quadratically with the amplitude of the applied aberration. Conversely, 
looking at points off-axis, the intensity is not constrained at a maximum value and so can vary 
linearly with the amplitude of the applied aberration. It is this linear sensitivity that results in 
significant distortion even at low levels of aberration. Of course the off-axis intensity of a planar 
beam is reduced relative to that on-axis and so any such distortions are not so significant in terms 
of the whole PSF. However, for singular beams the intensity in the PSF is concentrated off-axis 
and these distortions can become highly significant. We analyse these issues by calculating the 
sensitivity, Si, of a singular or normal beam to Zernike aberration mode, Z/, as the linear slope 
of the PSF intensity, /, with applied aberration amplitude, a{. 

The amplitude PSF, C/(p, t; ), is given by the Fourier transform of the unit radius pupil func- 
tion P(r,G): 
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where & [] is the Fourier transform over the unit circle from cylindrical co-ordinate space (r, 6) 
to cylindrical co-ordinate space (p, given as 
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The PSF intensity is then written as 
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Hence according to our definition the sensitivity Sj can be written as 
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Substituting equation eqQ]into eq|4]one obtains: 
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(4) 



(5) 



Zernike mode polynomials have the general form Z/(r, 6) = R n (r)cos(m6) or Zj(r, 6) = 
R n (r) sm(m6) with radial variation R n (r) and integer value m depending on the mode index, 
i. Note that circularly symmetric modes such as defocus and spherical aberration are repre- 
sented by modes of the first form with m = 0. Using the Fourier transform properties of such 
functions in polar co-ordinates fPflh we can write 
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(6) 



where J m is the Bessel function of the first kind of order m. Eq(5]can be simplified using eq|6j 
and basic trigonometry to give expressions for Si as follows. 
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With Z,(r, 0) = R„(r) cos(m0) 

st = (-l)?F/((-l)'G, 
Si = (_1)^F ; ((-1)'G, 

And similarly with Z/(r, 6) = R n (r) sin(m0) 
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Substituting m = in eqlTJ it can be seen that Si vanishes for all values of /, hence all modes 
are insensitive to spherical aberrations. For / = beams (no singularity) it can be seen that 
equations eq|7] and eq|9] yield Si = for even m. For odd m and / = beams, the sensitivity 
on-axis (p = 0) where the intensity in the PSF is highest, must also be zero from eq[l2l Hence 
a non-singular beam is only sensitive to odd aberrations, and then only off-axis where the 
intensity is low. Conversely for singular beams (/ ^ 0) and non- spherical aberrations (m ^ 0), 
Si always shows a cos(m^) or sin(m<^) functional form off-axis (and in some cases on-axis) 
where the power in the singular beam is concentrated. It is this result that demonstrates the 
linear sensitivity of PSF intensity to azimuthal aberrations for singular beams. 



3. Design of the sensor detector 



Discussion in the above section suggests that by analysing the azimuthal harmonic content of 
the PSF it should be possible to determine the azimuthal content of the aberrations present 
on the beam, which could then be used to provide a feedback signal for correction of those 
aberrations. Rather than using a whole image transform approach to achieve this, we suggest a 
segmented detector approach that would be more applicable to simple or low light level appli- 
cations. This form of detector would be easy to implement as a simple segmented photodiode, 
and with relatively few steps required for computation of the aberration signals, would be ap- 
plicable to real-time aberration correction applications. 
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Fig. 2. Gray scale images representing the magnitude of the sensitivity matrix, \Sij\ for 
(a) as singular / = 1 beam and (b) a non- singular / = beam. The displayed gray scale is 
normalised such that the maximum value, S55 with / = 1, is set to white, (c) The segmented 
sensor detector pattern. 



The sensor comprises twelve detectors, d n , each equal to a segment of angular width | radian 
of a reference circle as shown in fig Etc) and generating a signal 8 n in = 1 — > 12). The reference 
circle and its radius are termed as the sensor circle and sensor radius respectively. The sensor 
circle is symmetrically positioned in the focal plane of the aberrated singular beam. The detector 
signals from the sensor are split into two groups for each Zernike mode, namely a positive 
group, ?t + and a negative group, n~, as given in tableQ] 

AlJ ~ H8 n (13) 



dAij 



(14) 



The sensor signal, A^-, is defined as the response of the /th sensor, designed to detect Zernike 
mode Zj, to the presence of Zernike mode Zy. This is obtained by subtracting the total number 
of photons collected in the negative half from that collected in the positive half, and then nor- 
malizing by the total number of photons collected, as shown in eq{l3l (note that for Z\$ the n + 
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Table 1 . Geometry and performance parameters for the proposed wave-front sensor 



signals are counted twice). With this the sensor is able to detect 8 non-rotationally symmetric 
Zernike modes (Z/ = 5 6,7,8,9,io,i4,^J 15), referred to henceforth as the 8 sensor modes. We also 
define the sensitivity of the /th sensor signal to the jth mode, Zy, as the slope of the sensor 
signal vs strength of aberration, ay, at aj = 0, as shown in eqfFfl 

4. Numerical simulation of the sensor performance 

The performance of the sensor was simulated numerically using Fourier techniques according 
to eq[T]to determine an appropriate sensor circle radius. The sensitivity of the sensor for each 
sensor mode was computed as a function of sensor radius, normalised by the radius of the first 
zero of the Airy disk, p . The resulting self sensitivities, S#, are plotted as a function of the 
normalised sensor radius, p/p , in fig|2a). It can be seen that for different sensor modes the 
sensitivity saturates at a radius p = 3p . The self sensitivity, Su, at this radius is given in table[T] 
and we use this radius for all our subsequent calculations and experiments. The full sensitivity 
matrices were then calculated for the sensor against all Zernike modes up to j = 22 (second or- 
der spherical aberration) with both singular (/ = 1) and normal (/ = 0) beams shown graphically 
in fig J3a) and fig. Ob) respectively. The self sensitivities with singular beam are strong with a 
certain amount of cross sensitivity present, not least due to our chosen simplified method of de- 
tection and sensor signal calculation. As suggested by eq|7]and eqj9j the self sensitivities with 
normal beam vanish for sensor modes with an even value of m, and are typically significantly 
lower than those for the singular beam. 

The sensor signal was also calculated for each sensor mode over the range — 1 < a\ < 1 
and is shown in figJ3b). The sensor gives a variable performance depending on the particular 
mode being sensed, the best being for Z5 and (astigmatism). It is also noted that the effective 
linear range of the sensor is relatively short with saturation appearing between 0.5 < \at\ < 1.0. 
The range of the sensor is limited, but this is a result of and should be considered along side 
the high sensitivity produced. This high sensitivity, although over a short range, is particularly 
useful when used in closed loop control systems. 

Along with the sensor signal the associated error in the sensed aberration is estimated ana- 
lytically based on the shot noise in the sensed photon stream, assuming a mean number of N 
photons and a Gaussian approximation to the Poissonian noise statistics. This estimated error, 
Gu, is also shown in table [T] 

5. Experimental implementation 

A proof of principle system was constructed as shown in figlU Beam expanded and collimated 
light from a laser system operating at 532nm illuminates a reflective ferro-electric liquid crystal 




spatial light modulator (FLCSLM - CRL Opto SXGA). A computer generated hologram, de- 
signed using the technique described in OS [16]], is displayed on the device over a circular area 
of diameter 512 pixels. The hologram can be dynamically updated from a LAB VIEW program 
running on a PC and, along with lens L3 and iris I, effectively allows the accurate generation of 
arbitrarily shaped wavefronts across the 512 pixel system pupil. Light reflected from the FLC- 
SLM is directed towards lens L3 where the PSF of the system is realised in its back focal plane 
on the +1 diffracted order of the computer generated hologram. This beam is then imaged on 
to a fluorescence slide (FS), emitting in the red, with a 10X microscope objective (MO). A lens 
L4 forms an incoherent image of the resulting intensity PSF on the CCD camera (Hammamatsu 
ORCA ER) panel. A red filter (F) in front of the camera ensures that it receives only fluorescent 
light from the slide. This system allows the capture of clean images without coherent artefacts 
caused by scatter or reflections in the CCD camera. FigUa) shows a comparison between the 
experimentally generated and simulated PSFs for a singular (/ = +1) mode with 0.447Z; of 
various Zernike aberrations applied, clearly demonstrating the effectiveness of the system. 

In the system as realised the radius of the Airy disk was approximately p = 20 pixels and so 
a sensor radius of 60 pixels was chosen. Captured images were read by a LAB VIEW program 




Experiment (b) 




Iteration (t) 

Fig. 5. (a)Experimental (top) and simulated (bottom) intensity PSFs for the / = 1 singular 
beam aberrated with 0.447X Z^^^c,, 10,14,15 (left to right), (b) (92 KB) / = 1 PSFs (left) 
and corresponding normal / = PSFs (right) as a sequence (t=0— >4) during the correction 
procedure for a well corrected beam (t=0), aberrated by a glass slide, (c) (320 KB) Se- 
quence (t=0— >8) of / = 1 (top left), / = (top right) PSFs and the corresponding holograms 
(bottom) during the correction of a heavily aberrated beam, (d) Peak intensity of the / = 
PSFs as a Strehl ratio vs iteration, for the images shown in (b) and (c). 



running in the PC and the sensor signal was calculated by the program according to eq[T3l A 
small error is introduced in the sensed aberrations when matching the detector segments to the 
pixellated CCD array. With 1 1 ,096 pixels in the sensor circle we estimate this error to be less 
than 0.0022 radians for each of the sensor modes. Alignment of the detector with respect to the 
centre of the unaberrated beam is also a crucial requirement for the detection of correct signals. 
The lateral offset of the PSF can be considered in terms of the particular Zernike modes Z2 
and Z3 which represent tilts of an otherwise planar wavefront in the pupil plane of the lens. 
However, when using an / = beam figHb) shows that the sensor signals A7 and Ag become 
predominantly sensitive to Z2 and Z3 (S72 = S^3 = —0.92) and can therefore be used to centre 
the beam on the sensor circle. 

The performance of the sensor was evaluated experimentally by determining the sensor sig- 
nal for each sensor mode as the amplitude of that mode was varied. The results are plotted in 
fig 13c) showing clear agreement with the simulation results in fig Ob). 

Closed-loop control of aberrations was also demonstrated by using the detected sensor sig- 
nals as the feedback correction factor in a simple proportional closed loop control system. The 
movie in fig Ob) shows an initially well corrected beam at iteration t = that is aberrated at 
t = 1 by placing a glass microscope slide in the beam path. For this small, mainly astigmatic 
aberration the correction is made in a single step and in the following 3 iterations the PSF re- 
mains stable with only a residual drop in Strehl attributed to reflections from the glass slide. 
This Strehl ratio is plotted at time t = and subsequent times in figEfd). 

As a second demonstration of closed loop operation a random mix of the sensor modes, 
totalling (j)RMS = 1-4 radians, was applied to the system as an internal initial offset. The sub- 
sequent evolution of the system PSF along with the corresponding holograms is shown in the 



movie in figJ3c), and the Strehl ratio is plotted in figOd). Again the correction of the applied 
aberrations is clearly demonstrated in about 4 iterations, with the system exhibiting stability 
over the remaining 3 iterations. It is noted that in both these cases we considered only the diag- 
onal elements of the sensitivity matrix and made no effort to optimise the control loop for 
speed or stability. 

6. Conclusion 

We have shown that singular beams are ultra sensitive to non-rotationally symmetric aberra- 
tions, such that even well corrected optics produces significant distortion of the beam. We also 
show that this distortion manifests itself as an intensity modulation in the PSF of the system 
with the same azimuthal dependance as an applied Zernike aberration mode. We have suggested 
how this can form the basis of an aberration sensor that can measure those aberrations and ex- 
perimental results demonstrate the efficacy of our aberration sensor to correct non-rotationally 
symmetric aberrations. The aberration sensor also shows advantages in comparison to other 
sensing techniques in that there may be no need to include extra components in the system than 
are already there. Indeed it could even be considered to provide enough sufficient visual feed- 
back for the manual correction of aberrations that have been introduced by simple component 
misalignment. 
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